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In this study, we have found a new random ordered phase in isotropic models with many- 
body interactions. Spin correlations between neighboring planes are rigorously shown to 
form a long-range order, namely coplanar order, using a unitary transformation, and 
the phase transition of this new order has been analyzed on the bases of the mean-field 
theory and correlation identities. In the systems with regular 4-body interactions, the 
transition temperature T c is obtained as T c = (z — 2)J/k-g, and the field conjugate 
to this new order parameter is found to be H 2 . In contrast, the corresponding physical 
quantities in the systems with random 4-body interactions are given by T c = \/ z — 2J /k-g, 
and if 4 , respectively. Scaling forms of order parameters for regular or random 4-body 
interactions are expressed by the same scaling functions in the systems with regular or 
random 2-body interactions, respectively. Furthermore, we have obtained the nonlinear 
susceptibilities in the regular and random systems, where the coefficient Xnl of H 3 in 
the magnetization shows positive divergence in the regular model, while the coefficient 
X7 of if 7 in the magnetization shows negative divergence in the random model. 

Keywords: Many-body interaction; coplanar correlations; random ordered phase 



1. Introduction 

Recently, investigations of systems with many-body interactions are attractive even 
in the field of statistical mechanics of information, i.e., the learning theories us- 
ing higher-order Boltzmann machines^. Besides, some insula tors have been known 
as materials composed mainly of 4-body interactions^^], as shown in Eq.([T]). 
Study of phase transitions of systems with many-body interactions may be ex- 
pected to give some useful hints for solving the above information theoretic prob- 
lems. Wegner^l concluded that spin systems only with 4-body interactions have no 
spin ordered phase on the basis of the dual transformation method in 1971. How- 
ever, one of the present authors (M.S.) solved rigorously this 4-body spin model 
with extremel y an isotropic interactions ( J 2 = 0) in three dimensions using the a-r 
transformation^^. By this study of the extremely anisotropic model (we may call 
it no-ceiling model (fuki-nuke model in Japanese)) , Suzuki^ found the following 
interesting results, that is, (i) there exist an infinite number of degeneracies in the 
ground states, (ii) no spontaneous magnetization appears for all temperatures, (iii) 
spin correlations between neighboring planes (namely coplanar spin correlations) 
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form a long-range order below the transition point T c , (iv) the specific heat di- 
verges at T = T r - 

Savvidy et.all^2UQJ studied recently the same anisotropic 4-body spin model 
from a new view point of information theory and solved this mode l in dependently 
from Suzuki 6 using the same cr-r transformation. Castelnovo et.alP^ studied also 
the same models from an information theoretic interest on glass transition, and they 
conjectured that similar orders may appear in the isotropic models. It is not yet 
understood how the phase transition occurs and what configurations appear in the 
ordered phase of this isotropic model. 

In order to analyze the new random ordered phase, we consider here the following 
model 

% = — ^ JijkiSiSjSkSi — fisH^^^Si, (1) 

plaqucttes * 

only with 4-body interactions {Jijki}- In section 2, spin correlations between neigh- 
boring planes are rigorously shown to form a long-range order using the unitary 
transformation, and the phase transitions are analyzed on the basis of the mean- 
field theory and correlation identities. In the systems with regular 4-body inter- 
actions, the transition temperature T c is given by T c = (z — 2) J/fce- In contrast, 
it is given by T c = \fz — 2 J /fee in the systems with random 4-body interactions. 
Scaling forms of order parameters for regular or random 4-body interactions are 
expressed by the same scaling functions in the systems with regular or random 2- 
body interactions, respectively, as given in section 3. In section 4, we derived the 
nonlinear susceptibilities, where the coefficient Xni of H 3 has a positive divergence 
in the regular model, while the coefficient X7 of H 7 has a negative divergence in the 
random model. In conclusion, it is understood how the coplanar spin correlations 
form a random ordered phase, and how the phase transitions occur in systems with 
many-body interactions. 



2. Degeneracy and coplanar orders 

As was mentioned in the previous section, the "no-ceiling model" with 4-body in- 
teraction was solved exactly by one of the present authors (M.S.) 6 and was found 
to show very interesting behaviors in the low temperature phase below the transi- 
tion point T Ci o (which agrees with that of the 2d square Ising model). Fortunately, 
this is now found to play a crucial role in studying our isotropic model ([TJ. First, 
the transition point T c of this model is easily concluded to be higher than T c ,0: us- 
ing Griffiths inequalities^!. Furthermore, the peculiar behaviors of the " no-ceiling 
model" are also proved to be preserved in our isotropic model by introducing an 
infinite number of unitary transformations, as shown below. 

In general, using the Pauli operators {cf^ is id ^} on the grid points 
(ii,i2, ■ ■ ■ , id) defined in a d-dimensional space including orthogonal coordinate axes 



September 22, 2010 0:48 WSPC/INSTRUCTION FILE paper "arrive 



New Random Ordered Phase in Isotropic Models with Many-body Interactions 3 

Xi, X2, • ■ ■ , and Xd, the Hamiltonian with 4-body interactions {J\ m ... id ] is given by 

oy \ ^ \ ^ rim z z 

Z^i 2-^i il>— > id ,ii,...,i m ,..., i d ) (ii,. ..,«!,..., i m +l,-,<d) 

(i,m) (ii,-- ,id) 

x<T (i 1 ,...,i 1 +l,...,i m ,...,i i ) fi 1 ,...,i,+l,...,i m +l,...,i i )) ( 2 ) 
where the natural numbers l,m, and k are within the range [l,d] and the natural 
numbers {ik} are within the range [1,L] for the system-size L. We introduce the 
following unitary transformations 

II ^,...,4,.. A). ( 3 ) 

(h,— ,ik,--- ,id) 

which inverses all spins in the plane Xk — i 2 ■ The Hamiltonian % is invariant under 
any unitary transformation Uk, that is, the Hamiltonian commutes all the unitary 
operators Uk] UkWA^ 1 = H. Consequently, each order can appear independently 
in each plane. Therefore, there appear an infinite number of ordered phases in the 
spatially isotropic model, similarly to the spatially anisotropic modeH Especially, 
the number of degeneracy of the ground states corresponds to the number of possible 
numbers of the unitary transformations. Clearly, there exist d x L combinations of 
Xk and if as give possible unitary transformations. Thus, the ground state has 
2 Ld (<C 2 L ) degeneracies. Hence, it is shown that there exist the ordered states 
shown in Fig. [T] in the region T < T c , which appear independently in each plane. 

Then, we prove that similar configurations appear in other directions. The order 
parameters r?f fe fe = (S( iu ... <ikt ... ^S^,... ,i k +i,... ,i d )) ifor any (ii,-- - ,ifc,--- ,id)j are 
reexpressed in the operator forms rff£ = (of^ ... ik ... j^o?^ ... i k+ i ... j d ))- The spin 
correlation 

= ( CT (ii,.",i^".,i fc ,.",t < ,) cr f<i,.",t?."-,ifc+l.-" ) << l )) ^ 

for ii = if is transformed into 

«f s W = ((-^,..,i?,..A,...iri)(-°(li."-,iP,"-A+i,",«)> =C ( 5 ) 
under the unitary transformation 

II <>&,■•■,<?,•■ ,<-)' W 

(ill'-' ,*!)••• .id) 

which inverses all spins in a plane = if for / 7^ k. The relation jy?* = r]^ k is proved 
trivially for xi = ii ^ if , as well. Thus, the relationship 

C = (7) 
is proved for any ii. On the other hand, it is understood that 77^ can take each 
value ±|i7^f I independently of rff k since the coplanar orders perpendicular to x\ 

axis = (<%!,... ,»,,... .i,,)^,- ,»,+!,..., *„)) (for any (ii,--- meet the 



conditions 

^ =^^r 1 - ((-(ii,"^,..,^)^!,.,^!,.,^))) - -if, (8) 
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Fig. 1. This picture expresses the relationship of spin-pairs in the ground state, where the same 
color means the same relation (i.e. (+,+)(-,-) or (+,-)(-!+))■ One can prove in any dimensions that 
similar configurations appear in other directions. 



and 



7] 



(9) 



All the coplanar order specified by appears at the same temperature (seen 

in sec. 3). Thus, similar configurations shown in FigJT]are proved to appear in other 
directions. 



3. Analysis of order parameters 

In this section, we analyze the phase transit ions of the isotropic systems with 4-body 
interactions using the correlation identities™^ 1 . In order to apply these identities 
to our model, the Hamiltonian (p} is decomposed into the following two parts 



ft — Hi 



u' 



(10) 



where the partial Hamiltonian %' does not include the nearest neighboring spins Si 
and Sj . Then the partial Hamiltonian Hij including the spins Si and Sj is given in 
the form 



Hi 



-J(A 3 S Z + B 3 Sj + C 2 S i S j ) - »BH(Si + Sj 



(11) 
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Fig. 2. Spins interacting Si, Sj directly under the 4-body interactions. The black filled sites mean 
the spins in a plaquette. Notation A3 and B3 express the summation of (z — 2) productions of 
three spins, while C2 expresses the summation of (z — 2) productions of two spins. 

Here, A3 includes all products of three spins interacting directly with the spin Si. 
Similarly, B 3 includes all products of three spins interacting directly with the spin 
Sj . The term C2 , which is the most important term to analyze the coplanar order 
parameter, includes (z — 2) products of two spins interacting directly with the spins 
Si and Sj (Fig. [5]), where z denotes the numbe r of nearest neighbor spins. 

According to the correlation identitie d^^ , the spin correlation (SiSj) is derived 
easily as follows 

(SiSj) = ({SiSj)^) 

tanh KC 2 + t&nh(KA 3 + h) tanh(KB 3 + h) 



1 + tanh KC 2 t&nh(KA 3 + h) tanh( KB 3 + h) , 
(z~2)K(S k S l )+h 2 , (12) 



using decoupling approximations and relationships (A 3 ) = (S3) = (A3B3) = 0. 
These decoupling approximations are equivalent of mean- field approximations. The 
parameters h and K are defined by h = ^lbH/UbT and K = J/UbT, respectively. 
The self-consistency r\ = (SiSj) — (S^Si), which requests the translation symmetry 
of coplanar spin-correlations, yields the order parameter 77 for small H and for 
T > T c as follows 

v = ( s ^ = -& kB T- H ( z -2)r (13) 

Therefore we obtain the transition temperature T c = (z — 2) J/fce and the external 
field H 2 conjugate to the order parameter r\. 

When the 4-body interactions {±J} are distributed randomly, the number of 
frustrated unit cells is just the same as that of non- frustrated unit cells (Fig|3j) . 
Then, the order parameter ( is defined using the double average [(SiSj) 2 ], where 
[■ ■ ■ ] expresses the random average for the distribution of {±J}. Similarly to the 
regular model, the order parameter £ is obtained as 

C°c ; (14) 



Then we also obtain the transition temperature T c = \fz — 2 J/fce and the external 
field H A conjugate to the order parameter^ for random 4-body interactions. 

The above results are listed in Table(U There is no ergodic path^ by which the 
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Fig. 3. Configurations of random 4-body interactions on a three dimensional unit cell. This figure 
shows all the possible configurations of random 4-body interactions ±J on the developed figure of 
a three-dimensional unit cell. All unit cells are categorized according to the frustration, that is, 
non- frustrated unit cells arc categorized into " Group 1" , while frustrated unit cells are categorized 
into "Group 2". Same configurations under the inversion of all {Jijkl} are omitted in this figure. 



system can move from one ground state to another ground state only by the local 
correlation inversion in the systems with regular 4-body interactions, because regu- 
lar models have no frustration in any unit cell. For this reason, spin-correlations in 
the regular models need larger fluctuation for disordered states than in the random 
models. Therefore the transition temperature of regular systems is higher than that 
of random systems. 



Table 1. The table of order parameters, critical temperatures, conjugate external fields and univer- 
sality classes for four kinds of interactions. The critical temperatures of the random models are lower 
than those of regular models. The powers of conjugate external fields correspond to the number of 
spins included in the order parameters. The regular and the random models have their distinctive 
universality classes. 



Interactions order parameter k^T c field universality class(7, 1//3, 8) 

regular 2-body(Ferr ojE] m = (Si) zj H (1,2,3) 

random 2-body(SGj3EHII q = [(Si) 2 ] y^J H 2 (1,1,2) 

regular 4-body v = (SiSj) (z - 2) J H 2 (1,2,3) 

random 4-body C = [{SiSj) 2 ] -Jz~^2J H 4 (1,1,2) 
metal-insulator transitiorffSl (1 5 1,2) 
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Furthermore, scaling forms of the order parameters are also obtained near the 
transition point (T ^ T c ), by taking into account the nonlinear terms of equations 
of state derived by the correlation identities. The equation of state of the order 
parameter rj in systems with regular 4-body interactions is obtained as rf + 3trj — 
3h 2 = for t = (T — T c )/T c . Then the scaling form of the order parameter rj is 
derived in the form 

^ /2 /(^), (15) 

where the scaling function f(x) is obtained as the solution of the cubic equation 
/ + 3/ — 3a; = 0. This function f(x) is known as the scaling function of the 
ferromagnetic model constructed of regular 2-body interactions. In addition to that, 
the power 2 of h 2 in the scaling parameter corresponds to the number of spins 
included in the order parameter r\ = (SiSj). In a similar way, the equation of state 
of the order parameter £, namely ( 2 + tC, — 4r = 0, leads to the scaling form of C as 
follows 

£ = t f(^)'> /(*) = 5^/1 + 2^-1)- (16) 

Here, the scaling function f{x) is the same as that of the spin-glass model con- 
structed of random 2-body interactions^!. By the way, the power 4 of ft. 4 in the 
scaling parameter corresponds to the number of spins included in the order param- 
eterC= [<5,S,-) 2 ]. 

Additionally, the critical exponents (3, 7, 5 are obtained immediately from the 
equations of state (or the scaling forms of order parameters). Then the universality 
classes are expressed by the sets of the critical exponents (7, 1/(3,6). In the present 
case, the universality class of our regular 4-body model is given by (7,1//?, S) — 
(1, 2, 3), while that of random 4-body model is given by (7, 1/(3, 5) = (1, 1, 2). These 
universality classes are also listed in Tabl^TJ It is clearly understood that the uni- 
versality class of each model depends on whether spin interactions are regular or 
random. Especially, the universality class of the random 4-body model is the same as 
that of the spin-glass model (i.e., the random 2-body model) which was obtained by 
one of the present authors (M.S.) in The universality class of metal-insulator 

transitions belongs to this category, as was shown phenomenologically by March, 
Suzuki and Parrinellcff^. 



4. Nonlinear susceptibility 

In this section, the singularity of the nonlinear magnetic susceptibility^ is analyzed 
near the transition point T c , because it is useful in observing the coplanar random 
ordered states experimentally. For the case of regular 4-body interactions, the partial 
Hamiltonian Hi including the spin Si is obtained as 



Hi = -JA 3 S l - iL B HSi, 



(17) 
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using A 3 (A 3 is also used in cq. fTTj) ). According to the mean-field decoupling ap- 
proximation 

(SjSkSt) = (Sj)(S k Si) + (S^iSjSi) + (Si)(SjS k ) = 2m V (18) 

due to the fact that one of (SjSk),(SkSi) and (SjSi) is vanishing, we derive the 
following equation of state 

m = 6(z- 2)Kmri + h (19) 

from the correlation identity ($,-) = (t&nh(K A3 + h)} under the first-order ap- 
proximation with respect to m, 77 and h. Since the relations m ~ XoH and 
r\ r~j H 2 /(T — T c ) hold near the transition point, the magnetization m is obtained 
as follows 

m = X oH + xmH 3 + ■■■ ■ xm<x+jr^y, (20) 

which shows a positive divergence in the coefficient of H 3 . 

Similarly to the regular model, the nonlinear magnetic susceptibility in the ran- 
dom model is obtained as follows: 

m ^ xoH + XsH 3 + X5 H 5 + XrH 7 + ■ ■ ■ ; X 7 oc ~jT^Y' (21) 

which shows a negative divergence in the coefficient of H 7 , while X3 an d X5 are 
finite at T = T c . 

5. Summary and Conclusion 

We have found a random ordered phase in the isotropic (even regular) Ising mod- 
els with fo ur-b ody interactions, as was conjectured by Castelnovo, Chamon and 
Sherrington^^. Spin correlations between neighboring planes are rigorously shown 
to form a long-range order (coplanar order) using an infinite number of unitary 
transformations, and the detailed properties of this phase transition have been an- 
alyzed on the basis of the mean-field theory and correlation identities. 

Our main results are the following: (1) The order parameters of our models both 
for regular and random 4-body interactions are represented by the coplanar spin 
correlations. (2) The external field conjugate to the order parameter rj = (SiSj) is 
given by H 2 in the regular system. The external field conjugate to the order param- 
eter £ = [(SiSj) 2 ] is obtained as H A in the random system. (3) The scaling functions 
of the order parameters are the same as for 2-body interactions. (4) We have clar- 
ified the singularity of nonlinear susceptibilities in the regular or random systems 
by studying the magnetization in each system, as in ReflT] The coefficient of the 
H 3 form in the nonlinear susceptibilities shows a positive divergence in the regular 
model, and the coefficient of the term H 7 , for the first, shows a negative divergence 
in the random model. Finally we should remark that the nonlinear susceptibilities 
are crucially useful to study phase transitions experimentally or numerically on the 
systems with many body interactions. 
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Fig. Al. Coplanar correlation of spin correlations. The spins Si,Sj,Sk and Si included in the 
correlation function C(r) = ((SiSj)(Sf c Si)) are located as shown. 
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Appendix A. Correlation function 

Correlat ion fu nctions between coplanar spins are derived from the correlation 
identitiesPGl. The correlation function C(r) is defined by 

C{r) = ((SiS^SkSt)}, (A.l) 

where the spins Si,Sj,Sk and Si are located as shown in Fig lAll The correlation 
function C(r) for the Hamiltonian (fTT]) for H — is derived as follows 

z-2 

C(r) = (SkSiiSiSj)^) = (S k Si tanhXC 2 ) ~ X^{C(f + ae t ) + C{r- ae,)}, 

i=l 

(A.2) 

where a denotes the lattice constant and e*j denotes the unit vector parallel to the 
li-axis. In the second order-approximation of C(r ± aei) with respect to the lattice 
constant a, we obtain the differential equation 

x 7 i— 1 1 

which has the solution 

C(r) ~ ^ (A.4) 

for r = \r\. This solution C(r) shows the Ornstcin-Zcrnike type correlation function 
because of the mean-field approximations. Furthermore, the correlation length £ is 
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given by 

i= [Ka 2 /(l-(z-2)K)} 1/2 7 (A.5) 
as is easily seen from Eq. (|A.3[) . 
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